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Spatial  and  temporal  resolution  of  transition  to  turbulence  inside  the  free-shear  layers  of  two  yawed  cir¬ 
cular  cylinders  is  the  subject  of  the  present  investigation.  These  physics  were  resolved  using  the  large- 
eddy  simulation  (LES)  methodology.  An  0-type  grid  was  implemented  such  that  the  spatial  scales  of 
the  LES  computation  fully  resolved  the  energy  range  physics  of  the  shear  layers  at  Reynolds  number 
Reo  =  8000  based  on  the  cylinder  diameter.  The  two  test  cases  modeled  the  cylinder  span  skewed  at 
angles  45°  and  60°  from  the  horizontal  axis.  Observations  revealed  the  same  transition  process  as  the  nor¬ 
mal  cross-flow  state.  Soon  after  separation,  Tollmien-Schlichting  disturbances  were  predicted  that 
evolved  into  Kelvin-Helmholtz  (K-H)  eddies  before  absorption  by  the  large-scale  Karman-type  vortices. 
These  eddies  defaulted  to  a  spanwise  wavy  pattern  similar  to  a  normal  cross-flow  due  to  their  three- 
dimensional  instability.  No  mixed  modes  were  found  between  the  K-H  (Bloor)  and  Strouhal  frequencies. 
The  effect  of  yaw  angle  shortened  the  transition  process.  As  a  result,  peak  turbulence  levels  inside  the 
wake  formation  zone  approach  the  downstream  cylinder  periphery.  In  addition,  the  dimensionless  fre¬ 
quencies  of  the  K-H  eddies  lie  above  the  normal  cross-flow  relationship  as  formulated  by  Bloor 
(1964).  Disparity  between  the  yawed  and  normal  cross-flow  states  was  further  emphasized  by  the 
shear-layer  transition  characteristics.  Although  each  property  displayed  the  expected  exponential  growth 
during  transition  to  turbulence,  their  dimensionless  form  was  miss-aligned  with  those  of  the  normal 
cross-flow  case.  Based  on  the  present  evidence,  additional  simulations  (and/or  experimental  measure¬ 
ments)  are  necessary  to  form  conclusive  arguments  regarding  the  expected  behavior  of  the  transition 
characteristics  within  the  free-shear  layers  of  yawed  circular  cylinders. 

Published  by  Elsevier  Inc. 


1.  Introduction 

Understanding  transition  is  paramount  to  estimating  the 
streamwise  formation  of  large-scale  vortical  structures  in  the  near 
wake  of  bluff  bodies  such  as  circular  cylinders.  Under  normal 
cross-flow  conditions,  transition  to  turbulence  within  the  sepa¬ 
rated  shear  layers  is  dependent  essentially  on  Reynolds  number 
where  critical  defines  layers  truncated  to  the  turbulent  separation 
point.  Below  critical,  the  transition  process  involves  early  forma¬ 
tion  of  Tollmien-Schlichting  (TS)  waves  followed  by  secondary  ed¬ 
dies  that  ultimately  mature  into  the  much  larger  Karman-type 
structures  (see  Fig.  la).  This  progression  is  accelerated  when  the 
cylinder  is  oriented  at  angles  of  incidence  (yaw)  to  the  flow  free- 
stream  direction.  Consequently,  the  respective  time  scales  of  vor¬ 
tex  formation  are  concurrently  reduced  where  zero  is  reached 
once  the  cylinder  aligns  itself  with  the  freestream.  From  a  Naval 
systems  standpoint  involving  towed  communication  cables,  arrays, 
etc.,  yaw  causes  a  rise  in  certain  flow  characteristics  such  as  fluctu¬ 
ating  wall  pressure  that  leads  to  undesirable  self  and  radiated 
acoustics.  Moreover,  the  cylindrical  systems  risk  operational  fail- 
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ure  or  structural  fatigue  due  to  the  immediate  increase  in  oscillat¬ 
ing  lift  and  drag  forces. 

Unlike  transition,  the  effect  of  cylinder  yaw  on  shedding  fre¬ 
quency  is  well-understood  via  the  many  experimental  and  numer¬ 
ical  investigations  over  the  past  five  decades.  The  earliest  works 
(Hanson,  1966;  King,  1977;  Smith  et  al.,  1973;  Van  Atta,  1968)  fo¬ 
cused  on  verifying  the  independence  principle  where  the  axial  flow 
characteristics  must  be  decoupled  from  the  concurrent  cross-flow 
components  to  adequately  scale  the  vortex  shedding  frequency 
and  base  pressure  by  the  normal  velocity.  For  the  most  part,  this 
principle  fails  at  angles  beyond  10°  from  the  cylinder  normal  axis 
(Snarski  and  Jordan,  2001).  Contrary  to  the  formation  frequency, 
a  definable  argument  that  estimates  the  vortex  shedding  angle  as 
a  function  of  cylinder  yaw  is  still  at  large.  While  recent  studies 
indicate  a  sizeable  reduction  of  the  shedding  angle  at  large  cylinder 
yaw,  others  dispute  any  disparity  between  the  two. 

One  of  the  first  comprehensive  explorations  of  the  transition 
process  from  laminar  separation  to  turbulent  vortex  formation 
comes  from  Bloor  (1964).  Within  the  cylinder  separated  shear  lay¬ 
ers,  she  hypothesized  that  transition  initiates  with  the  appearance 
of  two-dimensional  small-scale  instabilities  whose  streamwise 
growth  amplifies  into  TS  waves.  These  waves  undergo  measureable 
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(a)  Transition  Process 


(b)  Secondary  (K-H)  Eddies 


Fig.  1.  Transition  process  to  vortex  formation  of  normal  flow  past  a  circular  cylinder  at  sub-critical  Reynolds  numbers. 


spanwise  distortions  that  accelerate  the  transition  process.  Bloor 
coined  the  early  disturbances  as  “transition  waves”  where  her 
hot-wire  anemometer  measurements  found  their  streamwise 
appearance  (after  separation)  to  be  inversely  proportional  to  the 
cylinder  Reynolds  number  (Reo);  Rcd  =  DUqIv  where  D  is  the  cylin¬ 
der  diameter,  Uq  is  the  freestream  velocity  and  v  is  the  kinematic 
viscosity.  This  detection  agreed  with  the  earlier  measurements 
by  Schiller  and  Linke  (1933)  where  they  noted  a  drop  in  the  mean 
transition  point  Xtr/D  =  1 .4-0.7  for  Rod  =  3500-8500.  Bloor  also 
noted  unique  dependence  of  the  transition-wave  frequency  (Jb) 
on  Rcd  when  scaled  by  the  formation  frequency  (Js)  of  the  large- 
scale  Karman  vortices.  She  formulated  the  relationship  (Jblfs)  oc 
Re^  given  /^  =  1/2  for  5000  <  Rod  <  25000.  This  1/2  power-law  was 
defended  by  Bloor  under  the  assumption  that  fb  scales  with  the 
boundary  layer  thickness  at  separation.  Later,  Unal  and  Rockwell 
(1988)  extended  the  lower-bound  detection  of  the  small-scale  dis¬ 
turbances  to  Rod  ^  1900  where  their  spectral  amplitudes  varied 
significantly  along  the  cylinder  span. 

Wei  and  Smith  (1986)  enhanced  our  physical  insight  of  the 
transition  process  to  turbulent  vortex  formation  for  circular  cylin¬ 
ders.  Like  Bloor,  they  detected  early  formation  of  small-scale  TS 
waves.  But  their  visualization  experiments  also  revealed  the  crea¬ 
tion  of  Kelvin-Helmholtz  (K-H)  spanwise  eddies  whose  shear- 
layer  evolution  originates  from  the  same  instability  mechanism 
described  by  Bloor.  Notably,  these  eddies  have  little  effect  on  the 
downstream  base  pressure  (Cpe)  inside  the  formation  zone.  Begin¬ 
ning  with  the  boundary  layer  at  separation,  Wei  and  Smith  re¬ 
ported  streamwise  transfer  of  the  respective  vorticity  to  the  K-H 
eddies  that  ultimately  was  absorbed  by  the  large-scale  Karman- 
type  structures.  However,  this  transfer  was  inconsistent  in  the 
spanwise  direction  which  lead  to  an  eventual  distortion  of  the  K- 
H  eddies  into  periodic  cellular  structures.  These  structures  aligned 
themselves  in  the  streamwise  direction  due  a  reverse  feedback 
mechanism  (downstream  to  upstream)  through  the  shear  layers. 
Wei  and  Smith  further  observed  breakdown  of  the  cellular  appear¬ 
ance  into  highly  irregular  structures  once  the  Reynolds  number 
reached  the  intermediate  range  (see  Fig.  lb). 

Besides  flow  visualization,  Wei  and  Smith  (1986)  conducted 
hot-wire  anemometer  measurements  where  they  reported  expo¬ 
nents  p  =  0.87  and  p  =  0.773,  respectively,  for  the  Bloor  power- 
law  expression.  They  reconciled  their  discrepancy  with  Bloor’s  Vi- 
law  by  scaling  fb  with  the  characteristic  momentum  thickness 
along  the  transition  length  rather  than  the  boundary  layer  thick¬ 
ness  at  laminar  separation.  Moreover,  they  revisited  Bloor’s  anal¬ 
ogy  over  a  wider  range  of  Reynolds  numbers  that  gave  the 
exponent  p  =  0.73.  However,  Wei  and  Smith  disbelieved  the  low¬ 
er-valued  exponent  based  on  the  premise  that  the  frequency  fb  is 
improperly  measured  when  using  anemometry.  Prasad  and  Wil¬ 
liamson  (1996)  compiled  their  anemometer  measurements  along 


with  those  of  Bloor  as  well  as  Wei  and  Smith  to  obtain  the  expo¬ 
nent  p  =  0.6742  over  cylinder  Reynolds  numbers  1.2  x  10^  ^  Rep  ^ 
4.5  X  10^.  Their  extensive  least-square  fit  provided  a  convincing 
argument  for  defining  the  scaled  transition  frequency  by  the  Bloor 
power-law  expression. 

Streamwise  pairing  of  the  K-H  eddies  was  not  visualized  by  Wei 
and  Smith,  nor  Unal  and  Rockwell  (1988),  in  the  cylinder  free-shear 
layers.  Consequently,  the  transition  frequency  fb  does  not  house 
sub-harmonics.  This  characteristic  allowed  Unal  and  Rockwell  to 
easily  quantify  dimensionless  flow  quantities  inside  the  shear  layer 
segment  that  owns  transition  to  turbulence.  For  example,  they 
found  an  exponential  streamwise  growth  of  the  dimensionless 
maximum  and  edge  velocity  fluctuations  as  well  as  the  integrated 
kinetic  energy  (see  Fig.  2).  Moreover,  each  quantity  collapsed  onto 
a  single  streamwise  variation  once  scaled  by  their  extrapolated  val¬ 
ues  at  boundary  layer  separation.  Their  growth  reaches  a  plateau 
once  transition  completes  which  disappears  at  the  critical  Rey¬ 
nolds  number. 

Few  numerical  studies  have  advanced  our  physical  understand¬ 
ing  of  the  cylinder  shear-layer  transition  process.  Braza  et  al. 
(1990)  computed  mixed  modes  of  nfs±mfb  similar  to  the  experi¬ 
mental  detections  of  Kourta  et  al.  (1987),  but  their  excessive  inte¬ 
ger  numbers  (n  and  m)  of  discrete  peaks  were  directly  attributed  to 
their  lack  of  adequately  resolving  the  spanwise  physics.  Further¬ 
more,  they  reported  streamwise  coalescence  of  the  K-H  eddies 
and  sub-harmonics  that  clearly  disagreed  with  the  experimental 
evidence.  Jordan  (2002)  investigated  the  transition  process  to  tur¬ 
bulence  at  Re^  =  8000  using  the  large-eddy  simulation  methodol¬ 
ogy.  By  ensuring  sufficient  spanwise  spatial  resolution,  he 
properly  predicted  the  Bloor  transition-wave  frequency  fb  includ¬ 
ing  two  mixed  modes  in  partial  agreement  with  the  spectral  mea¬ 
surements  of  Kourta  et  al.  His  computations  underscored  non¬ 
pairing  of  the  K-H  eddies  and  extended  the  streamwise  variations 
of  Unal  and  Rockwell  to  a  higher  Reynolds  number. 

At  present,  no  comprehensive  numerical  study  has  been  re¬ 
ported  about  the  transition  process  to  turbulence  for  the  yawed 
circular  cylinder  in  view  of  the  above  free-shear-layer  physics. 
Herein,  we  aim  to  begin  correcting  that  deficiency.  The  present  pa¬ 
per  examines  the  structural  and  quantitative  properties  of  transi¬ 
tion  inside  the  separated  shear  layer  at  Reo  =  8000  and  yaw 
angles  a  =  45°  and  a  =  60°  as  measured  from  the  horizontal  axis. 
Like  Jordan  (2002),  the  solution  methodology  is  the  large-eddy 
simulation  (LES).  The  fundamental  premise  of  an  LES  computation 
is  full  spatial  and  temporal  resolution  of  the  energy-containing 
scales  of  the  turbulence  motion.  The  finer  scales  (beneath  the  grid 
resolution)  are  commonly  universal  in  character  and  represented 
by  an  appropriate  sub-grid  scale  (SGS)  model.  Eor  the  cylinder 
free-shear  layers,  the  grid-scales  must  sufficiently  capture  early 
generation  of  the  small-scale  TS  waves  and  their  subsequent  local 
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streamwise  evolution  throughout  the  entire  transition  cycle.  Once 
resolved,  one  can  ascertain  direct  correlation  between  the  shear- 
layer  physics  of  the  yawed  cylinder  to  those  understood  for  the 
normal  cross-flow.  One  should  realize  however  that  this  study  only 
partially  addresses  the  overwhelming  range  of  possible  flow  phys¬ 
ics  that  manifest  within  the  separated  shear  layers  over  all  yaw  an¬ 
gles  and  Reynolds  numbers  below  critical. 


2.  Resolved  field  equations 


As  previously  noted,  the  LES  physics  comprise  the  resolved 
scales  by  the  local  grid  spacing  and  un-resolved  contributions  pro¬ 
vided  by  the  SGS  model.  Resolving  generation  and  growth  of  the 
shear-layer  TS  instabilities,  K-H  waves  and  fine-scale  physics  of 
the  Karman  vortex  formation  requires  grid  clustering  toward  the 
cylinder  surface.  This  stretched  structured  grid  necessitates  trans¬ 
forming  the  Cartesian  form  of  the  Navier-Stokes  equations  into  a 
curvilinear  coordinate  framework  =  c^,  rj,  z)  followed  by  a  math¬ 
ematical  grid  filter  to  acquire  the  final  LES  solution  set.  According 
to  Jordan  (1999),  this  spatial  order  of  coordinate  transformation, 
then  mathematical  filtering,  preserves  the  commutative  property 
between  the  differentiation  and  filter  operations.  The  former  trans¬ 
formation  operation  derives  a  set  of  direct  numerical  simulation 
(DNS)  equations  for  solution  in  the  computational  domain  while 
the  subsequent  filter  operation  over  uniform  curvilinear  spacing 
derives  the  respective  LES  formulation.  Eor  the  yawed  cylinder 
solutions,  we  will  assume  that  the  filter  width  is  synonymous  with 
local  grid  spacing.  The  final  grid-filtered  equations  in  curvilinear 
coordinates  become 


r-  n 

Continuity :  =  0 


(la) 


^  9\/lUj  dU'‘Ui 
Momentum  + 

di'‘  Redin'^ ^  oe. 


(lb) 


where  the  two  spatial  operations  redefine  the  real  SGS  stress  (erf)  in 
terms  of  the  resolved  Cartesian  (0,  v,w)  and  contravariant  (D,  V,  W) 
velocity  components;  g\  =  -  U'^Ui.  Each  contravariant  velocity 

component  is  evaluated  in  terms  of  their  resolved  counterparts 
using  the  definition  Uj.  While  the  overbar  in  this  defini¬ 

tion  denotes  the  filter  variable,  the  overtilde  symbolizes  implicit  fil¬ 
tering  of  the  metric  coefficients  (^f )  and  Jacobian  (v^)  via  their 
numerical  approximation. 


Temporal  evolution  of  the  transition  process  to  turbulence  was 
modeled  according  to  the  fractional  step  method.  The  convective 
terms  in  (1 )  were  spatially  approximated  by  a  sixth-order  compact 
stencil  in  the  circumferential  and  spanwise-periodic  z  directions 
with  an  optimized  fifth-order  template  in  the  wall-normal  vj  direc¬ 
tion  (see  Jordan,  2007).  These  terms  were  time-advanced  by  the 
second-order  explicit  Adams-Bashforth  method.  All  viscous  terms 
(including  the  SGS  model)  were  approximated  by  second-order 
central  differences  that  were  time-advance  by  the  implicit 
Crank-Nicolson  scheme  to  eliminate  the  associated  stability 
restriction  near  the  no-slip  walls.  Additional  information  regarding 
this  solution  methodology  was  presented  in  detail  by  Jordan 
(1999). 

3.  Dynamic  sub-grid-scale  model 

The  turbulent  scales  removed  by  the  spatial  filter  operation  are 
termed  the  sub-grid-scales  of  the  fluid  motion.  Eor  practical  appli¬ 
cations,  these  very  fine-scales  of  turbulence  are  most  commonly 
represented  by  the  dynamic  form  of  the  Smagorinsky  eddy  viscos¬ 
ity  relationship  (Smagorinsky,  1963;  Lilly,  1992).  This  form  will  de¬ 
fault  to  negligible  values  within  the  laminar  flow  regions  or 
turbulence  regions  of  high  grid  resolution.  Eurthermore,  dynamic 
evaluation  of  the  SGS  model  is  most  important  near  the  cylinder 
wall  segments  that  bound  the  upstream  side  of  the  vortex  forma¬ 
tion  region.  Toward  these  segments,  the  turbulent  stresses  typi¬ 
cally  display  an  asymptotic  character.  After  performing  the  two 
spatial  operations  described  above  for  the  grid-filtered  equations, 
the  dynamic  Smagorinsky  model  for  application  in  the  computa¬ 
tional  space  becomes 

af-l/3gfT,,  =  2Cd2|S|Sf  (2) 

where  the  single  dynamic  coefficient  C  is  evaluated  based  on  the  re¬ 
solved  scales.  The  metric  gf  is  a  grid-filtered  term  defined  as 
gf  =  The  turbulent  eddy  viscosity  (vt)  is  defined  by  the 

expression 

Vr  =  CZ2|S|  (3) 

where  |S|  =  and  A  is  the  grid-filter  width.  The  resolvable 

strain-rate  field  (Sf)  is  expressed  as  Sf  =  Sy  . 

After  following  the  least-squares  minimization  procedure  by 
Lilly  (1992),  dynamic  evaluation  of  the  SGS  model  coefficient  pro¬ 
ceeds  according  to  the  expression  (Jordan,  2001 ) 

C  =  ^  ,  (4) 

2A2M'' 
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where  the  tensor  Lf  is  the  contravariant  form  of  the  modified  Leon¬ 
ard  term  that  has  the  form 


The  Reynolds  stress  tensor  (if)  becomes  if  =  U%  -  U^Ui. 

The  tensor  Mf  is  a  modeled  contravariant  stress  that  is  defined 
in  terms  of  resolvable  tensors  as 

=  (6) 

wh^re  the  second  overbar  denotes  the  test-filter  width 

At]  At  =  2Ag  =  A  [0L  =  2). 

Unfortunately,  explicitly  filtering  the  local  instantaneous  strain- 
rates  by  the  test  filter  in  the  transformed  space  can  produce  nega¬ 
tive  values  for  the  model  coefficient.  These  contributions  were 
equally  seen  when  test  filtering  in  the  physical  domain  along 
non-uniform  spacing  (Jordan,  2001).  In  the  turbulent  energy 
spectrum,  positive  coefficients  denote  forward  scatter  while  the 
negative  values  symbolize  backscatter.  Inasmuch  as  energy  back- 
scatter  correlates  well  over  long  execution  times,  the  negative 
model  contributions  can  quickly  destabilize  the  LES  computation. 
Thus,  all  negative  coefficients  were  truncated  to  zero  in  the  present 
LES  computations  to  maintain  long-term  solution  stability. 

4.  Results  and  discussion 

In  the  following  section,  we  intend  to  investigate  the  funda¬ 
mental  physics  of  transition  to  turbulence  for  the  yawed  circular 
cylinder  at  test  angles  a  =  45°  and  a  =  60°  as  reference  to  the  hori¬ 
zontal  axis  (see  Eig.  3a).  More  specifically,  we  desire  to  discover  the 
initial  appearance  of  TS  waves  followed  by  the  streamwise  matu¬ 
rity  of  the  K-H  secondary  eddies.  These  physics  are  well-known 
to  form  within  the  free-shear  layers  immediately  after  laminar 
separation  under  normal  cross-flow  conditions.  Both  the  experi¬ 
mental  evidence  and  numerical  predictions  have  concurred  in 
identifying  this  process  ‘a-priori’  to  the  subsequent  formation  of 
the  downstream  turbulent  wake  at  Reynolds  numbers  below  crit¬ 
ical.  Notably,  this  particular  study  is  an  initial  numerical  attempt 
to  visualize  and  quantify  the  physical  impositions  on  the  shear- 
layer  transition  process  due  to  a  lateral  (spanwise)  velocity 
component. 

The  0-grid  generated  for  this  LES  application  housed 
241  X  141  X  64  points  in  the  peripheral  (s),  normal  (r),  and  span- 
wise  (z)  directions,  respectively.  This  grid  was  utilized  for  both 
LES  simulations.  Equidistant  clustering  was  invoked  along  the 
wall-normal  direction  towards  the  cylinder  surface  to  insure  that 


(Ar"”)  ^  2  for  the  first  field  point  throughout  the  wake  formation 
zone;  Ar'"  =  Aru^lv,  where  Ut;  is  the  averaged  wall  friction  velocity 
computed  along  the  cylinder  downstream  wall  within  the  forma¬ 
tion  zone.  The  wall-normal  outer  boundary  was  set  at 
^'max  =  ^  =  lOD.  Using  Kolmorgorov’s  microscale  (f/d)  as  experimen¬ 
tally  measured  by  Ong  and  Wallace  (1996)  and  the  relationship 
r]d  ~  Re“^^^,  we  can  ascertain  the  adequacy  of  grid  resolution  with¬ 
in  the  free-shear  layer  region.  This  resolution  is  illustrated  in 
Eig.  3b  where  the  profiles  depict  the  dimensionless  spatial  scales 
rjdjSx  and  rjdjSy  along  the  downstream  centerline;  and  dy  are 
the  centerline  grid  spacing.  Knowing  that  transition  completes 
within  x/D<1.5  at  this  Reynolds  number  (Jordan,  2002),  wave- 
number  cut-off  over  the  region  of  numerical  interest  is  essentially 
restricted  to  the  upper  bound  of  the  inertial  sub-range.  In  terms  of 
wall  units,  this  spacing  translates  to  A+  19  at  x/D?^  1.5  where 
A^  =  Asu-t/v.  Thus,  the  energy  scales  of  the  transition  process  were 
fully  resolved  by  the  present  computations.  In  the  spanwise  direc¬ 
tion,  a  uniform  fine-grid  spacing  of  Azf^O.OlOrr  (Aj  ?^21)  was 
implemented  to  assure  capture  of  the  spanwise  irregularities.  This 
spacing  resolves  the  spanwise  wavelength  (4)  of  the  secondary  ed¬ 
dies  by  at  least  five  points  (to  sixth-order  accuracy)  according  to 
the  experimental  measurements  by  Mansy  et  al.  (1994)  and  Wil¬ 
liamson  (1995);  4/D  -  20Re-^/2 

External  flow  conditions  involved  application  of  constant 
streamwise  and  spanwise  velocities  {Uoo  and  W^)  on  the  upstream 
side  of  the  outer  boundary  with  the  contravariant  form  of  Euler’s 
equation  (and  continuity  for  closure)  along  the  downstream  side; 
Uoo  =  Uosin(a)  and  =  Uocos(a).  The  Euler  boundary  condition 
suitably  exited  the  wake  structures  with  no  distortions  where  con¬ 
tinuity  specifically  enforced  the  incompressibility  constraint.  The 
resultant  of  the  constant  upstream  velocities  gave  Re^  =  8000  for 
both  yawed  cylinders  to  allow  direct  comparisons  to  the  normal 
cross-flow  simulations  by  Jordan  (2002)  at  the  same  Reynolds 
number;  Uo  =  (U^  +  W^).  Periodic  velocity  and  pressure  condi¬ 
tions  were  enforced  along  the  spanwise  boundaries  to  close  the 
kinematic  environment.  To  expedite  the  LES  computations,  both 
test  cases  were  initialized  by  the  instantaneous  variables  from  Jor¬ 
dan  (2002)  with  appropriate  adaptations  to  the  present  flow  condi¬ 
tions  and  grid  size.  A  low  time  step  (AT  =0.001)  was  fixed  to 
advance  the  LES  solutions  where  the  statistical  evidence  was  not 
collected  until  complete  elimination  of  all  flow  transients  as  given 
by  consistent  predictions  of  the  mean  lift  and  drag  forces. 

Before  exploring  the  free-shear-layer  physics  of  the  yawed  cyl¬ 
inder,  we  will  briefly  examine  the  respective  global  flow  properties 
of  the  a  =  45°  test  case  in  comparison  to  the  experimental  and 
computational  data  for  the  normal  cross-flow  state.  Table  1  lists 


(a)  Simulation  Conditions 


(b)  Spatial  Resolution 


Fig.  3.  LES  Conditions  of  Yaw,  spanwise  length  and  stretched  0-grid  spatial  resolution. 
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Table  1 

Comparison  list  of  several  global  parameters  to  the  experimental  data  of  normal  flow  past  a  circular  cylinder  (Jordan,  2002)  and  the  present  ^  =  45°  (^  =  90  -  a)  yaw  angle 
(Rod  =  8000). 


Parameter 

LES  ((/)  =  0°) 

Experiment  (0°) 

LES  ((/)  =  45°) 

Cl  (lift  coefficient) 

±0.115 

0.12  (Tadrist  et  al.,  1990) 

±0.119 

Cd  (drag  coefficient) 

1.06 

1.10  (White,  1974) 

1.04 

St  (Strouhal  number) 

0.204 

0.205  ±  0.002  (Norberg,  1993) 

0.207 

Up  (peak  velocity) 

0.46 

0.45  (Norberg,  1998) 

0.49 

CpB  (base  pressure) 

-1.08 

-1.06  (Williamson,  1995) 

-1.02 

fbifs  (Bloor  frequency) 

9.4 

10.1  ±1  (Prasad  and  Williamson,  1996) 

11.7 

6s  (separation  angle) 

92.2 

94  ±  2  (Son  and  Hanratty,  1969) 

93.2 

the  statistical  steady-state  properties  of  lift  coefficient  (Q),  the 
drag  coefficient  (Cd),  Strouhal  number  (SJ,  peak  streamwise  fluctu¬ 
ation  (Up),  formation  zone  base  pressure  (CpB),  scaled  Bloor  fre¬ 
quency  ifblfs)  and  separation  angle  {6s).  Several  subtle  differences 
are  clearly  evident  between  the  normal  and  yawed  cylinders.  For 
example,  the  coefficient  CpB  suggests  a  6%  reduction  in  base  pres¬ 
sure  with  little  reciprocating  affect  on  the  resultant  drag  coeffi¬ 
cient.  This  result  should  be  expected  in  view  of  the  surface 
pressure  response  to  a  non-zero  spanwise  velocity  component  in 
the  yawed  cylinder  case.  Conversely,  amplification  in  the  scaled 
Bloor  frequency /b//s  is  distinctly  apparent  in  Table  1  that  suggests 
earlier  formation  of  the  large-scale  Karman-type  vortices  coinci¬ 
dent  with  a  larger  spatial  and  temporal  decline  of  the  transition  cy¬ 
cle.  Evidence  supporting  reduction  of  the  transition  time  to 
turbulence  is  further  indicated  by  a  slight  increase  in  the  predicted 
lift  coefficient  Cl. 

While  the  mean  separation  angle  in  Table  1  appears  unaffected 
by  cylinder  yaw,  St  should  vary  according  to  the  independence 
principle  with  anticipated  variations.  The  LES  computations  reflect 
this  outcome  by  predicting  a  minor  increase  in  the  Strouhal  num¬ 
ber  for  the  a  =  45°  test  case.  But  a  better  representation  of  the 
shedding  behavior  is  illustrated  in  Eig.  4  where  the  Strouhal  num¬ 
bers  of  both  yawed  cylinders  are  scaled  by  the  normal  cross-fiow 
values.  This  figure  superimposes  the  present  results  as  well  as  pre¬ 
viously  reported  experimental  data  for  700  ^  Re^  ^1100  onto  the 
cos(a)  curve.  Like  the  experimental  evidence  at  much  lower  Re, 
both  LES  computations  gave  Strouhal  numbers  that  deviate  above 
the  independence  principle. 

Given  the  global  properties  of  the  present  test  cases,  we  can 
now  focus  our  attention  on  the  supportive  evidence  that  lie  inside 
the  separated  shear  layers.  One  of  the  first  concerns  of  interest  is 
the  evolution  and  spanwise  distortion  of  the  K-H  structures  them¬ 
selves  and  their  immediate  impact  on  the  near-term  formation  of 
the  Karman-type  features.  Jordan  (2002)  did  not  report  large-scale 


90-a 


Fig.  4.  Comparison  of  predicted  Strouhal  numbers  (a  =  45°  and  a  =  60°)  to  the 
experimental  evidence  and  independence  principle. 


spanwise  variability  of  the  secondary  eddies  for  the  normal  cross- 
flow  which  agreed  with  the  visualization  testing  of  Wei  and  Smith 
(1986)  at  their  higher  Reynolds  numbers  {Reo  >  3500).  This  obser¬ 
vation  is  only  partially  supported  by  the  present  a  =  60°  simula¬ 
tion.  As  illustrated  by  the  iso-surfaces  (2  = -2)  and  vorticity 
contours  in  Eig.  5,  the  K-H  eddies  are  well-organized  with  a  long 
wavelength  distortion  (2s  37rD/4)  that  is  aligned  in  the  stream- 

wise  direction.  The  parameter  2  as  used  to  identify  the  iso-surface 
structure,  is  the  second  eigenvalue  of  the  strain-rate  tensor  of  the 
instantaneous  turbulent  motion.  The  spanwise  wavy  pattern  of 
the  K-H  eddies  shown  in  Eig.  5  agrees  with  the  experimental  visu¬ 
alizations  reported  by  Lin  (1981).  At  comparable  Reynolds  num¬ 
bers  (based  on  the  normal  velocity  component),  Wei  and  Smith 
reported  these  small-scale  structures  to  possess  significant  irregu¬ 
larity  (more  turbulent)  with  intermittent  streamwise  alignment. 
One  can  also  see  in  Eig.  5  that  the  organized  K-H  structures  impact 
early  stages  of  vortex  roll-up.  While  these  eddies  are  visibly 
stretched  further  downstream,  their  spanwise  variability  is  mim¬ 
icked  by  the  initial  formation  of  the  larger-scale  Karman-type 
structures.  This  observation  agrees  with  our  basic  understanding 
of  the  latter  stages  of  transition  to  turbulence  where  the  Karman 
vortices  receive  a  majority  of  their  spanwise  vorticity  from  the 
K-H  eddies. 

A  snapshot  of  the  instantaneous  spanwise  vorticity  that  visibly 
shows  the  streamwise  evolution  of  the  cylinder  transition  process 
(Eig.  la)  within  the  upper  separated  shear  layer  is  illustrated  in 
Eig.  6  for  the  LES  test  case  a  =  45°.  Earliest  signs  of  the  two-dimen¬ 
sional  TS  instabilities  emerge  at  streamwise  position  x/D  0.26  as 
reference  to  the  separation  point.  Amplification  of  the  TS  waves 
into  K-H  eddies  is  plainly  visible  over  range  0.45  ^  x/D  ^  1.07 
where  the  upper  bound  marks  the  mid-point  of  a  fully  matured 
secondary  structure.  At  this  position,  this  eddy  is  speeding  through 
the  shear  layer  at  streamwise  velocity  u/Dq  0.56  or  ujlfoo  ~  0.80. 
A  second  K-H  eddy  is  observable  at  position  x/D  1.43  with  a 
higher  instantaneous  velocity  u/Dq  0.71.  Ostensibly,  this  velocity 
rise  is  rooted  in  the  adjacent  freestream  as  well  as  the  nearby  pull 
of  the  Karman  vortex  formation  which  together  manifests  a 
streamwise  increase  in  separation  gap  between  successive  K-H  ed¬ 
dies.  The  final  detectable  eddy  in  Eig.  6  owns  position  quantities  x/ 
D  1.58,  u/Do  0.21  and  v/Dq  0.23.  This  particular  eddy  has  al¬ 
ready  been  absorbed  by  the  much  larger-scale  Karman  vortex. 

A  second  set  of  snapshots  of  the  spanwise  vorticity  through  the 
formation  zone  and  bounding  shear  layers  are  shown  in  Eig.  7  for 
both  test  cases.  These  slices  depict  vorticity  contours  where  the 
values  in  the  parenthesis  denote  their  maximums,  minimums 
and  increments.  The  dimensionless  streamwise  and  spanwise 
velocities  (u,  w)  shown  over  each  K-H  vortex  represent  time-aver¬ 
aged  quantities  over  ten  transition  cycles  at  the  eddy  position 
(u  =  Umean/Lloo)-  According  to  these  averages,  the  K-H  eddies  accel¬ 
erate  in  both  the  streamwise  and  spanwise  directions  as  well  as 
mature  in  size;  du/dtieo  ~  0.12  and  du/dtUs  ^  0.32.  Their  stream- 
wise  growth  in  velocity  and  size  agrees  in  principle  with  the  exper- 
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(a)  Structure 


(b)  Vorticity  (dot-dash  box) 


Fig.  5.  Spanwise  structure  (A  =  -2)  and  vorticity  contours  (max.  10,  min.  0,  incr.  0.5)  for  yawed  circular  cylinder  at  Reo  =  8000  (a  =  60°,  ^  =  30°). 


Fig.  6.  Snapshot  of  the  spanwise  vorticity  contours  (max.-ll,  min.-16,  incr.  0.5 
and  max.  25,  min.  -25,  incr.  2.5)  showing  the  TS  instabilities  and  K-H  Eddies  within 
the  free-shear  layer  of  the  yaw  circular  cylinder  (oc  =  45°). 


imental  measurements  made  by  Unal  and  Rockwell  (1988).  In  par¬ 
ticular,  the  size  variation  is  attributed  to  their  exponential  increase 
in  turbulence  levels  with  downstream  distance.  A  new  contribu¬ 
tion  of  the  present  study  is  the  slight  downstream  acceleration  of 
the  K-H  structures  in  the  spanwise  direction;  dw/dtieo  ~  0.022 
and  dw/dt|45  0.017.  Although  these  levels  appear  to  correlate 
with  yaw  angle,  they  fall  well  below  their  relative  magnitudes  in 
the  streamwise  direction. 

Besides  their  spanwise  distortion,  the  streamwise  evolution  of 
the  K-H  structures  shown  in  Fig.  7  is  visibly  irregular.  We  note  that 
these  two  snapshots  of  the  instantaneous  shear-layer  physics  were 
chosen  specifically  to  reveal  this  inconsistency.  While  the  left  fig¬ 
ure  (a  =  60°)  suggests  an  integrated  time  scale  of  secondary  eddy 


generation  on  the  order  of  the  Strouhal  scaling,  the  right  snapshot 
(a  =  45°)  clearly  indicates  their  independence.  In  the  former  snap¬ 
shot,  one  can  easily  see  a  premature  Karman  vortex  that  has  al¬ 
ready  absorbed  smaller-scale  K-H  eddies  of  the  lower  shear- 
layer.  The  latter  snapshot  shows  K-H  eddies  (upper  shear  layer) 
stretched  and  displaced  rotationahy  as  they  near  the  larger-scale 
Karman-type  structure.  But  the  final  destiny  of  the  lower  eddy  in 
the  right  snapshot  is  somewhat  unclear  and  its  streamwise  evolu¬ 
tion  is  an  instance  of  strong  asymmetric  character.  One  should  note 
however  that  these  irregular  behaviors  are  sufficiently  infrequent 
such  that  the  upper  and  lower  transition  layers  share  the  same  sta¬ 
tistical  evidence. 

Another  observation  in  Fig.  7  is  the  streamwise  shortening  of 
secondary  eddy  spacing  with  increasing  yaw  angle;  Ax/DUs  = 
0.43  and  Ax/Djeo  =  0.40.  Evidence  supporting  this  characteristic  is 
the  probability  distributions  of  the  streamwise  K-H  eddy  gap 
(Ax/D)  while  passing  through  the  shear  layer  (see  Fig.  8).  Mean 
gap  lengths  are  =  0.47  and  ji  =  0.40  for  the  a  =  45°  and  a  =  60°  test 
cases,  respectively,  where  spacings  below  these  values  occur 
earlier  in  the  transition  phase.  This  gap  reduction  versus  yaw  angle 
should  be  expected  in  view  of  the  K-H  eddy  formation  frequency 
and  strength  responding  to  the  higher  wall  vorticity  at  laminar 
separation.  But  several  artifacts  of  yaw  angle  decline  are  earlier 
signs  of  turbulent  transition  and  vortex  roll-up,  shortening  of  the 
overall  formation  zone,  non-detection  of  K-H  eddy  coalescence. 


(a)  a  =  60^  (30,  -30, 2) 


(b)  a  =  45"  (32,  -28, 2) 


Fig.  7.  Snapshots  of  spanwise  vorticity  contours  (maximum,  minimum,  increment)  highlighting  the  secondary  (K-H)  eddies  for  yawed  circular  cylinder. 
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Fig.  8.  Probability  distributions  of  K-H  eddy  gap  distances  within  free-shear  layers 
of  two  yawed  cylinders;  F(Ax/D)  =  P(Ax/D)/P(Ax/D)max. 

and  closer  proximity  of  the  peak  velocity  fluctuations  to  the  down¬ 
stream  cylinder  periphery.  Evidence  of  the  latter  result  is  illus¬ 
trated  in  Fig.  9  for  the  streamwise  velocity  fluctuations  (<u'^>/ 
y^w  angle  a  =  45°  where  initial  detection  of  instabilities  oc¬ 
curs  soon  after  separation.  Interestingly,  this  turbulent  quantity  as 
well  as  its  streamwise  length  from  the  mean  separation  point 
(1.38D)  show  similar  values  as  the  normal  cross-flow  results  when 
scaled  by  the  normal/freestream  velocity  ratio  {UodUo).  Although 
not  shown,  we  note  that  this  simple  scaling  worked  equally  as  well 
for  the  a  =  60°  test  case. 

Knowing  the  downstream  limits  and  clear  appearance  of  K-H 
eddies  in  the  free-shear  layers  for  both  test  cases  involving  yawed 
cylinders,  we  can  now  survey  the  transition  process  to  turbulence 
from  Bloor’s  perspective.  First,  we  wish  to  distinguish  Bloor’s  tran¬ 
sition-wave  frequency  (Jb)  from  the  Strouhal  shedding  cycle  (Js)  fol¬ 
lowed  by  clear  detection  of  mixed  modes  on  the  order  nfs  ±  mfb. 
This  approach  is  illustrated  in  Fig.  10  by  plotting  the  downstream 
evolution  of  streamwise  energy  spectra  (Fuu)  for  the  a  =  45°  test 
case.  These  profiles  represent  spanwise  averages  over  six  Strouhal 
shedding  cycles  with  a  50%  overlap  (693  total  averages).  All  fre¬ 
quencies  and  energies  are  appropriately  scaled  by  the  cylinder 
diameter  (D)  and  the  resultant  freestream  velocity  (Ho).  The  line 
in  each  plot  denotes  the  5/3s  slope  of  the  inertial  sub-range.  Lastly, 
the  relative  downstream  positions  of  each  profile  are  given  by  the 
successive  shear  layer  cuts  in  Fig.  lOh. 

Emergence  of  the  TS  instabilities  was  found  in  the  energy  spec¬ 
tra  at  x/D  =  0.24  (Fig.  10a)  as  measured  from  the  mean  separation 
point.  This  location  occurs  about  25%  earlier  than  the  normal 
cross-flow  case  (x/D  =  0.32).  The  waves  themselves  hold  dimen¬ 
sionless  frequencies  near  1  Hz  where  the  spectral  peaks  below  this 
value  are  simply  harmonics  of/s.  Further  downstream  at  position  x/ 
D  =  0.33,  the  TS  instabilities  develop  into  K-H  vortices  where  loca¬ 
tion  x/D  =  0.44  marks  full  maturity.  The  dimensionless  Bloor  fre¬ 
quency  of  these  eddies  attains  the  distinct  value  /bD/Ho  =  1.71. 
But  further  downstream  between  regions  0.54  ^  x/D  ^  0.63,  this 
frequency  becomes  somewhat  broad  banded  which  is  contrary  to 


the  experimental  and  numerical  evidence  for  the  normal-cross- 
flow  state.  Moreover,  the  energy  peaks  are  comparatively  reduced 
within  this  range  and  the  expected  signs  of  mixed  modes  are 
clearly  absent.  This  latter  result  suggests  that  introduction  of  the 
velocity  component  (w),  which  essentially  acts  normal  to  the 
shear-layer  (x-y)  plane,  prohibits  measureable  interaction  be¬ 
tween  the  K-H  and  Strouhal  vortices.  The  turbulence  matures  rap¬ 
idly  within  this  range  that  substantially  disrupts  survivability  and 
coalescence  of  the  K-H  eddy  structure.  Still  further  downstream  at 
position  x/D  =  0.74  the  shear  layer  forms  an  inertial  sub-range 
where  the  small-scale  vortices  are  barely  discernible  in  the  respec¬ 
tive  velocity  spectra.  Transition  has  terminated  at  location  x/ 
D  =  0.95  where  the  shear  layer  became  fully  turbulent  according 
to  the  energy  spectra  in  Fig.  lOf.  Interestingly,  this  downstream 
limit  is  27%  shorter  than  the  normal  cross-flow  computation  (x/ 
D  =  1.31)  by  Jordan  (2002)  but  confirms  the  visualization  observa¬ 
tions  (x/D  1)  by  Szepessy  (1994). 

Recalling  the  apparent  disappearance  of  mixed  modes  within 
the  free-shear  layer  of  the  yawed  cylinder.  Fig.  11a  shows  a  sample 
energy  spectra  of  the  a  =  60°Computation  at  downstream  location 
x/D  =  0.53.  Again,  the  dimensionless  length  is  referenced  to  the 
mean  separation  point.  Like  the  previous  observations,  this  sample 
profile  indicates  that  the  spectral  peaks  below  the  Bloor  frequency 
{fbO/Uo  =  1.85)  are  harmonics  of  the  Strouhal  shedding  cycle  {fsDj 
Ho  =  0.179).  The  K-H  vortex  energy  levels  are  comparatively  low 
and  any  mixed  modes  were  plainly  absent  after  attempting  to  find 
suitable  integer  values  for  m  and  n.  As  illustrated  in  Fig.  11b,  this 
latter  result  may  explain  why  Bloor’s  frequency  scaling  does  not 
align  with  the  relationship  (fblfs)  oc  ReoP  for  both  test  cases.  The 
w-velocity  component  presumably  inhibits  consistent  coupling  be¬ 
tween  the  frequencies  fs  and  fb  that  disavows  a  simple  and  unique 
dependence  of  their  ratio  on  Reynolds  number.  However,  deviation 
from  this  frequency  scaling  appears  to  correlate  well  with  yaw  an¬ 
gle.  But  additional  evidence  would  be  necessary  to  form  a  conclu¬ 
sive  argument. 

As  a  final  look  at  the  shear-layer  physics  of  the  yawed  cylinder. 
Fig.  12  shows  the  downstream  distributions  of  dimensionless  inte¬ 
grated  kinetic  energy  {£')  and  momentum  thickness  (0/D)  at  four 
selected  spanwise  planes  (z/D)  for  the  a  =  45°  test  case;  namely, 
z/D  =  71/32,  71/4,  71/2  and  37r/4.  Both  distributions  are  evaluated 
across  the  shear  layer  and  represent  time  averages  over  55  cycles 
of  transition  to  turbulence.  Cylinder  yaw  was  recognized  in  the  en¬ 
ergy  distributions  by  modifying  the  original  relationship  of  Unal 
and  Rockwell  (1988); 

where  {6m)  is  the  mean  momentum  thickness  and  He  is  the  local 
edge  velocity  (see  Fig.  2).  With  reference  to  the  shear  layer,  6m  is 
taken  at  the  mid-point  between  the  computed  limits  of  transition. 
Like  the  normal  cross-flow  state,  these  shear  layer  properties 


(a)  90** 


(b)45“ 


Fig.  9.  Contours  of  scaled  streamwise  velocity  fluctuations  for  normal  (0.21,  0,  0.01 )  and  yawed  (0.24,  0,  0.024)  cylinders  (Reo  =  8000);  lengths  referenced  to  mean  separation 
point. 


496 


S A.  Jordan/ International  Journal  of  Heat  and  Fluid  Flow  31  (2010)  489-498 


display  an  exponential  evolution  throughout  the  transition  process.  wise  location.  Completion  of  the  transition  process  is  readily  visible 

This  evolution  is  clearly  consistent  but  not  independent  of  its  span-  by  the  abrupt  termination  of  this  streamwise  growth.  Both  profiles 


Fig.  10.  Scaled  spectral  energy  profiles  showing  streamwise  evolution  of  the  transition  process  through  the  cylinder  separated  shear  layer  (a  =  45°). 
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suggest  transition  to  turbulence  concludes  at  x/D  0.88  which 
coincides  with  the  energy  spectra  plotted  in  Fig.  lOg.  Similar  pro¬ 
files  of  the  a  =  60°  simulation  show  that  transition  ended  at  x/ 
D  -  0.92. 

Given  these  characteristic  distributions,  we  are  now  interested 
in  establishing  similarity  with  the  normal  cross-flow  case.  Follow¬ 
ing  the  earlier  work  of  Unal  and  Rockwell  (1988)  and  Jordan 
(2002),  both  properties  are  spanwise  averaged  over  the  streamwise 
limits  of  transition  with  the  resultant  distribution  of  £'  subse¬ 
quently  scaled  by  its  pseudo-intercept  at  x/D  =  0.  Results  of  this 


scaling  are  plotted  in  Fig.  13a  where  the  characteristic  streamwise 
length  over  the  transition  limits  was  changed  to  x/0m-  Only  a  cur¬ 
sory  look  is  necessary  to  conclude  that  the  yawed  cylinder  distri¬ 
butions  are  plainly  dissimilar  from  the  normal  cross-flow.  The 
relative  slopes  suggests  that  the  growth  of  E'  diminishes  with 
yaw  angle,  but  additional  test  cases  would  be  required  to  explore 
a  definable  association. 

On  the  contrary,  the  downstream  distributions  of  OjD  coincide 
for  both  yawed  cylinders  as  shown  in  Fig.  13b  with  similar  expo¬ 
nential  growth  rates  as  the  normal  cross-flow.  Although  intriguing. 


(a)  Spectral  Energy 


Fig.  11.  Profile  of  scaled  streamwise  kinetic  energy  at  x/D  =  0.53  for  yaw  cylinder  (a  =  60°)  and  scaled  Bloor  frequency  versus  Reynolds  number;  ^Bloor  (1964),  ^Wei  and 
Smith  (1986),  ^Prasad  and  Williamson  (1996). 


Fig.  12.  Streamwise  evolution  of  the  integrated  kinetic  energy  and  momentum  thickness  through  the  separated  shear.  Layer  at  four  spanwise  planes  of  a  yawed  cylinder 
(a  =  45°). 


Fig.  13.  Streamwise  evolution  of  the  integrated  kinetic  energy  and  momentum  thickness  (spanwise  averaged)  through  the  separated  shear  layer. 
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Fig.  14.  Mean  growth  rate  of  momentum  thickness  through  transition  of  normal 
and  yawed  cylinder  (Re  based  on  normal  velocity). 

this  discovery  demands  further  computations  (and/or  experiments) 
at  varying  conditions  to  attain  conclusive  proof.  This  statement  is 
clearly  supported  in  Fig.  14  where  the  mean  growth  rate  of  6  as 
evaluated  through  the  transition  region  is  plotted  versus  Reynolds 
number  (Re)  based  on  the  normal  freestream  velocity  (doo)-  Both 
rates  deviate  substantially  from  the  characteristic  curve  for  the 
normal  cross-flow  case  over  the  range  of  intermediate  Reynolds 
numbers. 

5.  Final  remarks 

The  present  LES  computations  of  transition  to  turbulence  within 
the  separated  shear  layers  of  two  yawed  circular  cylinders  reveal 
very  interesting  results  when  compared  to  each  other  as  well 
as  to  the  normal  cross-flow  state.  At  the  Reynolds  number 
ReD  =  8000  (below  critical)  and  yaw  angles  a  =  60°  and  a  =  45°  (as 
measured  from  the  horizontal  axis),  transition  to  turbulence  owns 
the  same  fundamental  characteristics  as  the  normal  cross-flow. 
The  transition  process  initiates  with  similar  appearances  of  Tollm- 
ien-Schlichting  (TS)  waves  soon  after  laminar  separation  that  rap¬ 
idly  matured  into  Kelvin-Helmholtz  (K-H)  vortices  before 
concluding  with  formation  of  the  large-scale  Karman  structures. 
Taken  collectively  with  the  previous  experimental  evidence,  one 
is  drawn  to  strongly  argue  for  this  transition  process  to  be  a  ubiq¬ 
uitous  mechanism  within  the  free-shear  layers  of  circular  cylinders 
or  bluff  bodies  in  general.  This  mechanism  seems  to  easily  survive 
cross-flows  that  act  normal  to  the  shear  layer.  However,  a  lower- 
bound  of  this  transition  process  must  exist  as  the  cylinder  span  ap¬ 
proaches  alignment  with  the  freestream  flow.  This  lower  limit 
would  distinguish  a  time  scale  for  the  Karman  vortex  formation 
where  the  K-H  eddies  would  simply  not  evolve  within  the  sepa¬ 
rated  shear  layer.  We  saw  this  omission  in  the  snapshot  of  the 
LES  computation  at  a  =  60°.  But  a  convincing  answer  to  this  ques¬ 
tion  would  require  extensive  supplementary  evidence. 

The  characteristic  physics  within  the  free-shear  layers  of  the 
two  yawed  cylinders  followed  the  same  downstream  progression 
as  the  normal  cross-flow  condition.  The  integrated  kinetic  energies 
and  momentum  thicknesses  grew  exponentially  downstream  until 
the  free-shear  layer  became  fully  turbulent.  Distinguishing  their 
spanwise  variability  from  the  normal  cross-flow  case  may  be  an 
excellent  subject  for  later  investigations.  But  by  resembling  the 
spanwise-averaged  normal  cross-flow  progression,  these  charac¬ 
teristics  are  excellent  omnipresent  indicators  for  deflning  the  spa¬ 
tial  limits  of  transition.  Although  not  shown,  the  fluctuating 


properties  inside  the  free-shear  layer  of  these  yawed  cylinders  dis¬ 
played  equivalent  trends.  But  their  streamwise  growth  rates  did 
not  align  with  the  normal  cross-flow  state  as  was  discovered  for 
the  integrated  kinetic  energy  and  momentum  thickness.  As  men¬ 
tioned  above,  further  data  is  necessary  to  establish  tangible  behav¬ 
iors  of  the  free-shear  layer  properties  for  the  yawed  cylinder. 
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